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Abstract 



Copositive and completely positive matrices play an increas- 
ingly important role in Applied Mathematics, namely as a key 
I concept for approximating NP-hard optimization problems. The 

' cone of copositive matrices of a given order and the cone of com- 

pletely positive matrices of the same order are dual to each other 
with respect to the standard scalar product on the space of sym- 
metric matrices. This paper establishes some new relations be- 
• tween orthogonal pairs of such matrices lying on the boundary 

of either cone. As a consequence, we can establish an improve- 
I ment on the upper bound of the cp-rank of completely positive 

■ matrices of general order, and a further improvement for such 

. matrices of order six. 
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1 Introduction 



In this article we consider completely positive matrices and their cp-rank, 
as well as copositive matrices. An n x n matrix M is said to be completely 
positive if there exists a nonnegative (not necessarily square) matrix V such 
that M = VV~^ . An n xn matrix A is said to be copositive if x''^Ax > for 
every nonnegative vector x S ffi!J: . The completely positive matrices of order 
n form a cone, C*, dual to the cone of copositive matrices of that order, Cn- 
Both cones are central in the rapidly evolving field of copositive optimization 
which links discrete and continuous optimization, and has numerous real- 
world applications. For recent surveys and structured bibliographies, we 
refer to [H [3 [H [11] , and for a fundamental text book to [3j. 

A main motivation for this paper was the study of cp-rank: A given 
completely positive matrix M always has many factorizations M = VV~^ , 
where y is a nonnegative matrix, and the cp-rank of M, cpr M, is the min- 
imum number of columns in such a nonnegative factor V (for completeness, 
we define cpr M = if M is a square zero matrix and cpr M = oo if M 
is not completely positive). Determining the maximum possible cp-rank of 
n X n completely positive matrices, 

Pn ■= max {cpr M : M is a completely positive nx n matrix} , 

is still an open problem for large n (up to now, for n > 6; only recently 
j55 = 6 has been established [IS])- It is known [21 Theorem 3.3] that 



For n G {2,3}, there exist simple proofs of ([T]), but already for n = 4, the 
argument is quite involved [3]. For n > 5, it is known that 



but whether the lower bound d„ is in fact equal to p„ is still unknown for 
large n. This is the famous Drew-Johnson-Loewy (DJL) conjecture [lOj . 
The above upper bound on p„ on the right-hand side follows, for example, 
from the so-called Barioli-Berman [1] bound: Let 

hr ■= max {cpr M : M is a completely positive matrix with rank M = r} , 
then for r > 3 



Some evidence in support of the DJL conjecture is found in (TUl [3 12 dS], 
see also [3l Section 3.3]. The DJL conjecture has recently been proved for 
n = 5 [18], but the cp-rank problem is still not fully resolved. Not only is it 



Pn = 



n if 77, < 4 . 



(1) 






(3) 
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not known whether the DJL conjecture holds, but the best upper bound on 
Pn for n > 6 remained, for over a decade, Two main results of this paper 
are a reduction of the upper bound on p„ in the bracket ([2]) for general n 
and a further reduction in case of n = 6. To obtain these results, we use 
[18\ Thm.3.4], which guarantees that p„ is attained (also) at a nonsingular 
matrix on the boundary of the cone of C*. We also complement this result 
here by studying, for every possible cp-rank 1 < k < Pn, where in C* is 
cp-rank k attained. 

Each matrix on the boundary of the cone C* is orthogonal to a matrix 
on the boundary of the cone C„ (in fact, to a matrix generating an extreme 
ray of that cone). Thus to improve the bound on p„ we consider pairs of 
matrices, M S C* and A £ Cn, that are orthogonal to each other in the 
standard scalar product of matrices. This leads also to some results that 
are not directly related to the cp-rank problem, and are of interest in their 
own right. 

The paper is organized as follows: after introducing basic concepts and 
terminology, we show, in Section [2l some important orthogonality and di- 
agonal dominance results. Section [3] is devoted to the study of extreme 
copositive matrices of low rank, while Section H] deals with genericity of the 
property of having a fixed cp-rank within the completely positive cone. Sec- 
tion [5] presents improvements of upper bounds on the cp-rank for matrices 
of general order, and a further tightening of this bound for order six is put 
forward in the final Section [6l 

Some notation and terminology: let ej be the ith column vector of the 
n X n identity matrix In- The nonnegative orthant is denoted by M" . For a 
vector X G M'^, the support of x is denoted by cr(x), i.e., 

(t(x) = {i : Xi > 0} . 

The set of nonnegative nxp matrices is denoted by M!|^^^. A matrix A € 
MJ]^^ is called positive if minAj,- > 0. (Note that a completely positive 

matrix is not necessarily positive, since it may have zero entries.) A matrix 
A G R'^^"' is diagonally dominant if \Aii\ > ^ \ Aij\ for every 1 < i < n, and 

is strictly diagonally dominant if all these n inequalities are strict. For two 
square matrices A, B we denote 



A®B 



A O 
B 



where O is a suitable (possibly rectangular) zero matrix. 

By Sn we denote the space of real symmetric n x n matrices, and 
by T^n tlie cone of symmetric psd matrices, Vn = {X G 5„, : X ^ 0}. 
The cone of nonnegative matrices in Sn is denoted by Mn-, i-e., Mn = 
Sn n M^^". The scalar product of two matrices U,V of same order is 
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{U,V) := trace (U^V) 



UV = [vi,...,vp] G 



Rl""", then 

the factorization M = VV~^ is equivalent to M = Ylf=i^i'^J ■ refer 
to this sum as a cp decomposition. When p = cpr M we say that the cp 
decomposition is minimal (the cp factorization is minimal). 

By K° we denote the relative interior of a convex set K, dK is the 
boundary of that set. For a convex cone ext K denotes the set of all 
elements in K who generate extreme rays of K. 

Both the copositive cone C„ and the completely positive cone C*, are 
pointed closed convex cones with nonempty interior. As mentioned above, 
the copositive cone Cfi and, in particular, its extremal rays, are important 
for the study of the cp-rank as any matrix on the boundary 5C* of C* 
is orthogonal to an extremal ray of C„,. However, characterization of the 
extremal rays of C„ for n > 5 is itself a major open problem in the study of 
Cn- The explicit characterization of extremal rays of C5 was completed by 
Hildebrand [T^ only recently, and this work was essential for the arguments 
in [18]. One extremal ray of C5 is generated by the so-called Horn matrix 



H 



1 1 

-1 1 

1 -1 

-1 1 

1 -1 



which historically was the first copositive matrix detected outside of Vn +J^n 
(here n = 5) [7]; attribution to Alfred Horn can be found in [13]. In the 
sequel, a matrix A is said to be in the orbit of a matrix B \l A = DP^ BPD, 
where D is a positive-definite diagonal matrix and P a permutation matrix. 
The Horn orbit consists of all matrices in the orbit of if; obviously, each 
matrix in the Horn orbit also generates an extremal ray of C5. Finally, we 
address any extremal matrix in C5 which is neither in the Horn orbit nor in 
U As as a Hildebrand matrix; see description in [181 Thm.4.2]. 



2 Orthogonality and diagonal dominance results 

In this section we consider copositive and completely positive matrices which 
are orthogonal to each other. The following theorem will be used in this 
section to point out a property of matrices on the boundary of the copositive 
cone, and also later to reduce the upper bound for p„. 

Theorem 2.1 Let M G C* be orthogonal to A ^ Cn, and let M = X]j'=i ^i^J 
be any cp decomposition of M . Then 

(a) For every 1 < i < n the i-th column of M is orthogonal to the i-th 
column of A. 



3 



(b) If i € o'(xj) for every 1 < j < p, then the i-th column of A is in the 
nuUspace of M. 

Proof. The scalar product of the i-th column of M and the i-th column 
of A is 

p 

ejMAei = ^ eJxjxjAei = ^ (^J^j) (^^J ^^i) ■ (4) 

j = l j-iea{xj) 

(a) li i = i, the right-hand side in (j3D is 

j:iGo-(xj) 

Since each Xj is in and satisfies xjAxj = we have [^Xj]^, = for all 
k G cr(xj) [18\ Rem. 3. 2, (3.1)]. In particular [j4xj],. = for all j's in the sum 
above. Hence ej MAei = 0. 

(b) Let z = Aei. Suppose i £ o"(xj) for all j G {1, • • • ,p}, then as above 
[Axj] . = for every j, and thus by dH 

p p 
[Mz], = ejMAe, = (e/x,) (xj^e,) = (e/x,) [Ax,]^ = 

for every ^ G {1, . . . , n}. □ 



Before we proceed, we note an interesting implication about coposi- 
tive matrices on the boundary dCn- It is well known, and obvious by a 
Gershgorin-type argument, that singular matrices (e.g. those on dVn) can- 
not be strictly diagonally dominant. For matrices on dCn we show that 
some form of "anti-diagonal dominance" can be established: 

Corollary 2.1 Let A G 5C„. If A ± M £ C* and M has no zero rows, then 
A is in the orbit of some A G dCn which satisfies 

An < E \Aij\ for all i . 

Proof. If ^ = O, this is trivial. Else, for A G 3C„ \ {O}, we select 
M G C* n A-^ \ {0} and choose Uj G M?J_ where p = cpr M, 1 < i < n, such 
that Mij = ujuj (in other words, M is the Gram matrix of ui, . . . ,Un). 
Next we select i so that ||uj|| > ||uj|| for all j. Then, since M ^ O, we have 
||uj|| > and 

^Aijujuj = {Aei)'^{Mei) - AuWuif = -AuWuif , 
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by Theorem l2.1l fa). Passing to absolute values and applying Cauchy-Schwarz 
as well as the triangle inequality, we get 



Scaling M by a positive-definite diagonal matrix D so that M = DAI D has 
diag M = e, we get that A = D~^AD~^ satisfies A _L M, so that the above 



According to a result by Kaykobad [15] , any symmetric diagonally dom- 
inant matrix in Nn is already completely positive. Triggered by above, we 
could ask whether indeed these matrices are in the interior of C* . The answer 
is negative, a certificate being I„: matrices in the interior of C* are neces- 
sarily positive. In general, being nonsingular and positive is not a sufficient 
condition for an n x n completely positive matrix to be in However, 
we can prove the following 

Theorem 2.2 Let M he an n y. n symmetric matrix, n > 3. If M is 
diagonally dominant and positive, then M G [C*]°. 

Proof. Let and e = [1, . . . , 1]^ G M", J„ = ee""" denote the all ones n x n 
matrix, and let // := miuj j- Mjj > 0. Consider M' := M — fiJn- Since 
n > 3, M' € Mn is strictly diagonally dominant and therefore completely 
positive and nonsingular. So we can put M = VV~^ where vi = y^e, and 
the remaining Vj come from the cp factorization of M'. By Dickinson's 
characterization [8] of [C*]°, the assertion is proved. □ 

Note that for n = 2 there exist diagonally dominant matrices that are 
not in e.g., J2, which is singular, and therefore on 862- 

3 Extreme copositive matrices of low rank 

A € ext Cn n A^n, then there is at most one positive entry on or above the 
diagonal. If this entry is on the diagonal, we have rank A = 1 and A £ Vn- 
If the positive entry is off the diagonal, then A is in the orbit of the matrix 
-£'12 = GieJ + e2ej , and hence is of rank two. Next we will sharpen these 
assertions, basically dropping the nonnegativity assumption on A. We will 
need the following auxiliary result, on the role of zero entries on the diagonal 
of an extreme copositive matrix: 

Lemma 3.1 Suppose that A € ext C„ \Mn can be decomposed as 




inequality holds for all 



i, if we replace A with A. 



□ 



A = 



S R 
Q 



with S £ Sk and diag Q = o G R' 



A; > 1. 
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Then R and Q are zero matrices (of suitable orders) and S G ext Ck \ {O}. 

Proof. Since diag Q = o and j4 G C„, we deduce Q £ m(^"-'=)>^(«-^) and 
R G M^^^""'^^ Further, since A ^ Mn, S G R''^'' has at least one negative 
element. Thus k > 1 and S ^ O. We conclude that 



A = 



' s 


o ■ 


+ 


o 


R ' 





o 




Q _ 



where the rightmost matrix has no negative entries and therefore is copos- 
itive. As S ^ O, extremality of A implies that both Q and R have to be 
zero matrices, and A = S (B O as well as extremality of S in follows. □ 

We can now prove: 
Theorem 3.1 Let A G ext C„,. Then 

(a) rank A = 1 if and only if A is positive-semidefinite. 

(b) rank A = 2 if and only if A is in the orbit of E12 = eiej + 6261. 

Proof. The if parts are obvious. For the only if: 

(a) If rank A = 1 then, since A is symmetric, A = zbxx"'' for some x G M". 
Since the diagonal entries of A are nonnegative, A = xx^. 

(b) Suppose rank ^4 = 2; if A G Mn, then the result follows directly. So 
suppose A has a negative entry. By extremality A ^ Vn- Hence A 
must be indefinite, i.e. of the form A = uu^ — vv^ = xy"*^ + yx""" 
(take, e.g., x = ^(u + v) and y = u — v). For any z G we have 

< z'^Az = 2(x"^z)(y'^z) , 

in particular Xiyi > for all i. Put r := {i : Xiyi ^ 0}. Then Xiyi > 
for all i £ T. By permuting rows and columns if necessary, we 
may assume that r = {!,..., A;}, and A can be decomposed as in 
Lemma |3.H yielding A = S (B O. So we may without loss of generality 
assume that Xiyi > holds for all i, by investigating 5 instead of A. 
By diagonal scaling we may now further assume that Xiyi = 1 for all i. 
Now, if both X, y G or both — x, — y G M" , we again would arrive 
at ^ G Mn- So we are done if we reduce the assumption Xj > > xj 
ad absurdum. To this end, consider the 2x2 block corresponding to 
these two indices {i,j}, putting t = ^ < 0: 



t + 



t + 
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Copositivity of this 2x2 matrix is equivalent to the condition t + j> 
—2, which upon muhiphcation with t < amounts to (t + 1)^ < or 
t = —1. So ah positive entries of x are equal, say a, and all negative 
entries of x equal —a. Hence y is a multiple of x and rank A = 1 < 2, 
a contradiction. 

□ 



Remark 3.1 Theorem \ 3. 1\ implies that for n < 4 each matrix in ext Cn\{0} 
has rank 1 or 2. The characterization of the extreme copositive 5x5 matrices 
implies that each matrix in ext C5 \ {O} has rank 1,2, or 5 Jj^[ /. What are 
the possible ranks of matrices in ext C„, n > 6? Note that if there exists 
A G ext Cn of rank k, then there exist also matrices in ext Cn+i of rank k 
(A(BO, to name one). 

We proceed with an immediate consequence for positive and nonsingular 
matrices M G 

Corollary 3.1 If M e dC*\{dVn U dMn), then any A G ext C„nA/^\{0} 
satisfies rank A > 3. Moreover, no principal suhmatrix S of A can he in the 
orbit of E12. 

Proof. Because M is assumed to be positive, we know A ^ Mn] similarly, 
since M is nonsingular, we conclude A ^ Vn- Therefore Theorem 13.11 implies 
rank ^ > 3. Next suppose a principal submatrix S ^ O oi A were in the 
orbit of £^12 and thus has diag S = o. Then A can be decomposed into 



A = 



' s 


■ 


+ 





R ' 










Q _ 



where R has no negative entries and Q is copositive. Hence the rightmost 
matrix is copositive, and (by extremality of A and S ^ O) therefore must 
be the zero matrix. It follows rank A = rank S ®0 = rank S = 2, but then 
Theorem 13.1( b) yields the contradiction A G dMn- CD 



4 Genericity of matrices of fixed cp rank and order 

We now turn to the study of the cp-rank of matrices in C*. In this section 
we consider the question whether or not having a certain cp rank is robust 
within all cp matrices of order n or not. 

First we observe that every possible cp-rank is attained at some matrix 
on the boundary: 
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Proposition 4.1 For every 1 < k < pn there exists a matrix Mk G dC* 
such that cpr = k. 

Proof. By [THl Thm.3.4] there exists a matrix M G dC* with cpr M = pn- 
Let M = X]j=i be a minimal cp-decomposition of M, and let = 

Ylj=i^j'^J ■ Then cpr < A;, and strict inequality is impossible, because 
it would contradict the minimality of the cp-decomposition of M. That is, 
cpr Mk = k. Since M is on the boundary of C*, there exists A E dCn \ {O} 
such that M is orthogonal to A. Then VjvJ _L A for every j, and thus 
Mk -LA, and therefore Mk € ^C*. □ 



However, it is interesting to find out whether there are also interior 
matrices having a prescribed cp-rank, and whether this property is robust. 
For this purpose, we denote the set of completely positive matrices of order 
n with cp-rank exactly equal to k by 

C^k = {M G C: : cpr M = k} . 

The extreme case /c = p„ is easy: as shown in [181 Cor. 2. 5], the set C* p^ 
contains an open set, and thus, 7^ 0. To prove this for all other 

k, we need a result which may also be of independent interest. Beforehand 
note that every M £ [C*^p„Y a factorization M = VV^ where ^ > 
has Pn columns, and by the Diir and Still characterization of [0^]° |12] it 
is easy to deduce that there exists a factorization M = WW^ where W is 
positive (and has rank n). However, this does not necessarily imply that 
there is a factorization M = VV~^ where V > has pn columns and all of 
these columns are positive. 

Proposition 4.2 There is always a matrix M G [C^,p„]° such that there 
exists a positive n x pn matrix V with M = VV~^ . 

Proof. Let Mq be some matrix in the interior of C* ^^. As in jl8j . let 
V S ° be its Perron- Frobenius eigenvector to the eigenvalue A > 0. If 
Mq = VoV(^ where Vq is a nonnegative n x p„-matrix, then no column Voe^ 
of Vq is zero. Therefore (V^^v)j = v^VoOj > for all i, in other words, 
i := Vf^v G [WlY. From Av = Mqv = Vo^V^v = Vqx it follows that 
X := x/A is a positive vector with v = Vqx. For small e > it follows from 
the choice of Mq that M = Mq -\- evv^ also has cpr (M) = pn. Now 

M = Mo + evv^ = VoVj + {Vo^)e{Vo^y = Vo{In + exx^)Fo^ . (5) 

For 5 = (Vl + exTx - l)/||x|| > 0, define C = In + (Jxx'^. Then = /„ 
exx""^ and V = VqC = Vq + 6{Vox)x^ is positive (as before, Vqx £ []R"]°). 
By dSD we obtain VV^ = VoC^Vf^ = M. □ 
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Theorem 4.1 

Proof. For k < n it follows from cpr (A/) > rank (M) that C* ^ is 
contained in the set of matrices with rank at most k and thus its interior 
is empty. We now show that [C* 7^ if n < A: < To this end, 
Proposition 14.21 ensures we can select a matrix M = VV~^ G with 
a positive n x pn matrix V = [vi, . . . , Vp^]. As M G [C*]° C P„°, we 
have rank V = n and without loss of generality, let the first n columns 
{vi, . . . , v„} of V be linearly independent. Now, let any k with n < k < pn 
be given and consider the matrix 

k 

i=i 

Obviously cpr M < k. On the other hand, cpr M < k would contra- 
dict the minimality of the factorization M = VV~^ , so cpr M = k. Let 
V := [vi,...,v„] and V := [v„+i, . . . , v^]. Then V is positive and non- 
singular square, so by |12^ Thm.2.3], we have M = VV + VV~^ G l^nT ■ 
Next consider the singular value decomposition of = UiT,U2 with suit- 
able orthonormal n x n matrices Ui and U2 and a positive-definite diag- 
onal n X n matrix S. Let l/JSnU^ be the set of all matrices of order n 
which result from premultiplying a symmetric matrix Z by Uj and post- 
multiplying it by Ui . Consider the map T : l/JSnU^ — )■ 5„, defined by 
J"(Ay) := (y + AV)(V + AV)'^. The derivative of J" at AF = is given 
by the Lyapunov operator 

Cy : UjSnU^ Sn with Cy{AV) = {AV)V^ + V{AVy . 

Given a symmetric right hand side R, solving Cy{Uj ZUj) = R for a sym- 
metric matrix Z is equivalent to 

UjZU^Ui^U2 + Ujm^UiZU2 = R 
ZT. + T.Z = Uj RU2. 

Evidently, this is uniquely solvable for a symmetric Z so that by the inverse 
function theorem, is invertible in an open neighborhood of AV = 0, and 
the inverse function satisfies V+AV > in this neighborhood, by continuity. 
Summarizing, for any (symmetric) matrix M in an open neighborhood of 
M there exists a positive n x n perturbation matrix 

V = V + F~^{M -VV^) 

of F, such that M = VV^ + VV^ = VV^ + Ei=n+i VjvJ G Q, which 
establishes cpr M < k. But we know from |18l Cor 2.5] that all matrices 
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M € C* which are sufficiently close to M have cpr M > k, so we conclude 
M G C* ^, hence M is an inner point of C* ^, and the results follow. □ 



5 New bounds for the cp-rank 

In this section we prove that the known upper bound bn on the cp-rank of 
n X n matrices can be reduced, for every n > 6. For n = 6 we reduce the 
bound further in the next section. First, we combine the idea of [17j with 
Theorem 12. II to show that p„ is strictly less than bn for every n > 3. 

Theorem 5.1 For n > 2, if A £ dCn has k > 2 positive diagonal elements, 
and M £ C* is orthogonal to A, then cpr M <bn — k + 1. 

Proof. We may assume that An > for i £ {1, . . . ,k}. Let 

£ = |S G cSn : ejBAei = , ah i G {1, . . . , . 

Then {vv'^ : v G M*^ and v'^Av = 0} C £ by Theorem [211 The subspace £ 
is isomorphic to the solution space of the homogenous system of k equations 
in variables bij, 1 < i < j < n, 

Aiibii + ^ ^ Aijbji + ^ ^ Aijbij = , i G {1, . . . , fc} . 
j<i i<j 

Since the diagonal matrix with An, i = 1, . . . ,k, on the diagonal is a subma- 
trix of the coefficients matrix, the rank of the coefficients matrix is k. Thus 
dim£ = ("+^) -k. Next suppose M G C* is orthogonal to A G dCn- Then 
M G conv { vv"*" : v G and v'''^v = 0} which is a convex cone contained 
in £, and by Caratheodory's theorem cpr M < dim£ = C^^^) — k. □ 

Thus for certain completely positive matrices on dC* we get the following 
bound on the cp-rank: 

Corollary 5.1 For n> 5, if A £ dCn \ Mn , and M £ C* is orthogonal to 
A, then cpr M < 6„ - 4. 

Proof. We may assume A £ ext Cn- If A is positive-semidefinite it 
follows from orthogonality and M ^ that rank (M) < n — 1 and thus, 
by ([3]), cpr (M) < < 6^ — 4. We now assume that A is indefinite. By 
Lemma 13.11 A has at least 5 positive diagonal entries, since otherwise we 
would have A = S (S O, where S* is a copositive matrix of order at most 4, 
so that S £ V4 + Mi, therefore A £ Vn + A/"™, and by extremality, A would 
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be either positive-semidefinite or nonnegative. 



□ 



It is beneficial to introduce a cp-rank bound for positive matrices at the 

boundary of C*: Let p* := max < cpr M : M G dC* , miuMj,- > >. 

I *J J 

Theorem 5.2 For n > 2, there exists a positive matrix M £ dC^ \ dMn 
such that 

Pn < cpr M + 1 . 

Hence we get 

P*n<Pn<P*n + l; (6) 

the right inequality is an equality for n G {2, 3} whereas the left inequality is 
an equality for n G {4, 5}. 

Proof. Let M £ [C*]° be a matrix such that cpr M = Pn [M Cor.2.5]. 
Let 5 > be such that M = M - (5e„eJ G Clearly, pn = cpr M < 

cpr M + 1 . Since M has positive off-diagonal entries in the last row and it 
is positive-semidefinite, we have Mnn > and thus M is positive, so that 
cpr M < p*. Hence p* < Pn < + 1- The last assertions follow from 
62 = 2 = p3 — 1 and from the fact that there are singular positive matrices 
M G dC* with cpr M = p„ for n G {4, 5} [El Rem.2.1,Cor.4.1]. □ 



Remark 5.1 For n < 4, the matrix M in Theorem \5.S\ is necessarily sin- 
gular. Thus for n = 2, rank M = 1 and cpr M = 1, which yields p2 < 2, 
a tight bound. If n £ {3,4} then rank M = n — 1, and thus cpr M < hn-i, 
and Pn < (2) ■ For n = 3 this gives ps < 3, which is also a tight bound. But 
for n = 4 this yields p4 < 6, which is not tight. For n = 5, it turns out 
that though p^ = is attained at a singular matrix fT^ Cor. 4-1], we have 
P5 = 6 < (2) = 10. Still, for n > 6 we get an improvement of the known 
bound Pn ^ bn, and for n = 6 we will improve it further below. 

By the above, for the first time we have a proof that bn is not a tight 
upper bound on the cp-rank of completely positive matrices of any order 
n > 5. More precisely: 

Corollary 5.2 For n > 5, we have 

Pn<bn- ^ and Pn<K-'i- 

Proof. Let M G 9C* \ dMn be a positive matrix with cpr M = p* . M 
is orthogonal to a matrix A G ext Cn- As M is positive, A Nn- By 
Corollary 15. H cpr M <bn — ^, and Pn < — 3 by Theorem 15. 2[ □ 
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Remark 5.2 For n = 6 this bound is pg < 17, but may be slightly improved. 
This case is studied in Section\^ 

Beforehand we note a further result vahd for arbitrary order. 

Theorem 5.3 If M G C* has a zero entry, i.e, Ms dJ\fn, then 

cpr M < 2p„_i . 

Proof. We may and do suppose that Mi„, = 0. Let M = X^f^;^ WjW^ 
be a cp-decomposition of M. Define := {i £ {!,..., : 1 G cr(wj)} 
and ^2 := {!,..., p} \ ^i, as well as Mj := X^jg^ . w-jW^^, for j = 1,2. 
Then Mi = M{ and M2 = M^^, where M[,M^ are matrices in C*_i. 
The result follows from cpr M = cpr (Mi + M2) < cpr Mi + cpr M2 = 
cpr M( + cpr Mg < 2p„_i. □ 



Remark 5.3 Note that for all n > 6, we have 2dn-i < 6n — 3, the bound 
from Corollary \5.2[ However, compared with the upper bound from Theo- 
rem \5.1l we have 2pn~i > 2d„_i > bn — k + 1, whenever k > ^, so that 
Theorem \5.3\ is interesting only for small k. 

6 Cp-rank of matrices of order six 

In this section we improve the upper bound on the cp-rank of completely 
positive matrices of order 6. First we consider matrices in dC^ which are 
orthogonal to S* 0, where S £ ext C5 is either in the orbit of the 5x5 
Horn matrix H or a Hildebrand matrix. Below, we use the symbol " ffl " to 
denote the sum modulo 5 of two elements in {1, 2, 3, 4, 5}. 

Proposition 6.1 Let S be either in the orbit of the 5x5 Horn matrix H 
or a Hildebrand matrix. Suppose that M £ Cq is orthogonal to S (BO. Then 
cpr M < 15. 

Proof. If M = VV~^ , V > 0, then each column of y is a nonnegative linear 
combination of three vectors, eg, e^ + effli and ejogi + 64^2 [IBl Thm.4.4]. 
Let 

W = [ei+ e2|e2 + esjes + e4|e4 + esjes + eijee] . 

Then V = WX where each column of X has support of at most 3 elements, 
contained in a set of the form {i, i ffl 1, 6}. For each 1 < i < 5, let Xi consist 
of the columns of X whose support is contained in {i, i ffl 1, 6}. Then, up to 
permutations of rows and columns, XiXj = 1^ with 1^ E C3 so that 

cpr XiXj = cpr < p3 = 3 . 
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5 5 

Therefore cpr XX'^ = cpr XiXj < cpr XiXj < 15. □ 

i=l i=l 



Theorem 6.1 pg < 15. 

Proof. By [H Thm.3.4], we know pQ = cpr M for some M £ dQ \ OVq. 
Moreover, if M had a zero entry, we get from Theorem 15.31 that cpr M < 
2p5 = 12. Suppose now that M e dC^\ {dVe U OMq). Then Corollary O 
gives rank A > 3 for all A S ext Cg H M-*- \ {O}, and at least one such A 
exists as M G OCq. Now either all diagonal elements of A are positive, in 
which case by Theorem 15.11 cpr M < — 5 = 15, or yl has at least one 
zero on the diagonal. By Lemma 13. H A = S ® O with S £ ext C5 . Since 
rank S = rank ^ > 3, we conclude that S is either in the orbit of iif or a 
Hildebrand matrix. Then Proposition 16.11 gives cpr M < 15, and the claim 
is proved. □ 

We thus cut the bracket for pQ in about half, since = 20 and = 9. 
The same argument could be used also for n S {7,8}, but it would not fur- 
ther improve upon the bounds yielded already by the general improvement 
in Corollarv l5.2[ 
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